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ABSTRACT 

In two phase flow systems line losses comprise frictional and mo- 
mentum pressure drops. For design purposes, it would be desirable to 
estimate the line losses employing a one-dimensional calculation. In 

• • f 

this paper two methods for computing one-dimensional momentum flux at a 
test section discharge station are compared to the experimental value for 
a range of two^-phase flow conditions. The one-dimensional homogeneous 
model appears to be more accurate generally in predicting the momentum 
than the variable slip model. 

INTRODUCTION 

The pressure drop in a two phase flow system comprises frictional 
and momentum pressure drop. For a two phase system in which the fluid 
is being heated, the latter can be the most significant because of 
rapidly changing momentum accompanying void fraction increases. In this 
paper we are going to examine the. validity of predicting momentum pres- 
sure drop by the use of a simple one-dimensional calculation. 

Many analytical and experimental investigations have been devoted 
to studying two-phase pressure drop in a one-dimensional system. The 
most widely used analyses are those presented in references 1 and 2. 
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The possibilities of so many types of flow patterns in which the gaseous 
and liquid phases flow at different relative velocities has made it im- 
possible to. separate the momentum and frictional pressure drops to one 
universal theory. Consequently, there have been a variety of analytical 
models suggested which have been evaluated for specific two phase flow 
regimes. The simplest of the flow models has been the homogeneous , in 
which it is assumed that both phases are uniformly distributed across 
the diameter of the passage and both flow at the same velocity (no slip) . 
The homogeneous model is usually treated as the reference one because of 
the relative ease with which it can be applied analytically. The more 
complicated models delineate possible variations in the void distribution 
and in the relative velocity between the phases (variable slip). While 
the more complicated models can relate more realistically to actual 
flows, there is always the question of how closely the simple homogeneous 
model will approximate certain properties of the flow - mainly pressure 
drop. Such an appraisal awaits more experimental information than has 
been generated thus far in which void, momentum pressure drop, total 
pressure drop, and flow rates are measured independently; 

In reference 3,' the authors measured the momentum flux discharge 
from a two component or a two phase flow in a pipe. For the two compo- 
nent case, they metered the air and water discharges; for the one- 
component case, they used a heat balance to estimate the liquid and 
vapor flow rates. The average void fraction was determined by isolating 
a section of transparent pipe by means of quick closing valves. One of 
the principle conclusions of reference 3 was that the one-dimensional 
modeling using measured void fractions (and variable slip) underpredicted 
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the system momentum and the homogeneous model, independent of the meas- 
ured void and with no slip, yields a better estimate of momentum. The 
results did raise questions about the homogeneous model applicability ! to. 
the data at the low quality regime. The authors raised the possibility 
of flow oscillations contributing to the disagreement between data and 
calculation. 

High velocity data pressure drop reported in reference 4 and dis- 
cussed in reference 5 indicate that one-dimensional modeling is quite 
acceptable at high system velocities. A composite of the experimental 
results of references 3 and 4 shows a definite velocity (Reynolds num- ••• 
ber) effect on the one-dimensional applicability. This effect is over 
and above that which is experienced in single phase flow as the flow 
transitions from a parabolic laminar velocity profile to a power law 
turbulent profile. This particular velocity dependence is a reflection 
of the system flow patterns and transitions (bubbly, slug, annular, mist, 

etc. ) . 

i 

Herein, one of the principle objects will be to perform a similar 
comparison between experimental data with the one-dimensional no slip 
(homogeneous) and variable slip representations for a range of velocities 
and void fractions.. 

Void fraction, flow rates of the gas and liquid species, and momen- 
tum discharge from the exit of the test section will all be measured in- 
dependently. The void' fraction measurement is made at a station just up- 
stream of. the test section discharge using the standard "one-shot" gamma 
attenuation technique. "One shot" means that the gamma beam is fixed at 
one position for all the tests. 
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The experimental momentum pressure drop results will be compared to 
the variable slip and the homogeneous model computations. The ranges of 
conditions wherein these models are suitable and wherein they are inac- 
curate will be discussed. The principle experimental information con- 
tained herein relates to the discharge momentum and the average void 
measured at the test section exit. 

APPARATUS 

The apparatus is diagramed in figure 1. The principal components 
comprise a test section, a nitrogen supply, and a water supply. The 
test section is a reinforced flattened tube with a gas-water mixer on 
one end for generating the cwo component flow. Flattening the test 
section confined the flow to an approximate one-dimensional pattern and 
increased the dimension in the direction of the gamma radiation, so as 
to increase the sensitivity of that measurement. The discharge from the 
open end of the test section was caught in a momentum cage enabling a 
measurement of the average force related to the momentum of the fluid 
discharge. A calibrated load or thrust cell was used to record the 
force. Figure 2 shows the construction of the momentum cage. • It is an 
assembly of aluminum plates arranged to turn the flow discharged from 
the test section through an angle of 90°. The cage is free-hanging, 
but is linked to the load cell behind it. 

As is shown in figure 1, just upstream of . the discharge point gamma 
radiation attenuation was used to measure the "one shot" void fraction. 
The gamma radiation was detected with a photomultiplier tube; the output 
was recorded on a strip-type recorder. A simple procedure in which a 
range of void blockages (lucite inserts with nearly the same attenuation 
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as water) were introduced into the gamma beam allowed the output of the- 
photomultiplier tube to be checked, A linear reproducible curve shown in 
figure 3 resulted. The abscissa of figure 3, which is the void fraction 
determined from gamma attenuation, was determined from the simple ex- 
pression 

In. $/<(>,. 

a = - — - — j— ( 1 ) 

meas In 

where <p Q , <j>^, and <p are the empty, full, and two-phase signals re- . . 
ceived from the photomultiplier sensor, respectively. Equation (1) was 
derived by Hooker and Popper (ref. 6). The derivation of equation (1) 
is based on the premise that the attenuation of a gamma beam is an expo- 
nential function of the absorption thickness. Henry (ref. 4) describes 
briefly other assumption inherent in equation (1) and estimates the 
errors associated with the use of this equation. Figure 3 shows good 
agreement between equation (1) and calibrated plastic voids which were 
introduced into the test section. During each run the all liquid and 
all gas conditions were compared. The signal strength ratios were con- 
sistently 1.29 which assured the operators that the void measuring sys- 
tem was on calibration. 

Before the water and the nitrogen were introduced into the mixer, 
their flow rates and temperatures were monitored separately. Chromel- 
alumel thermocouples were used for the gas and water temperature meas- 
urements. The orifice flowmeters which were used for the flow measure- 
ment were calibrated in a standards laboratory. The pressure drop 
across the orifices was sensed by strain-gage type pressure transducers. 
It is estimated that the flow rates were known to within ±2 percent. 
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In a two phase flow system there is always the question of relative 
stability of the flow. Strip chart records of all the flow transducer- ~ 
measurements were made, and the mean values were used to compute the 
flow rate. In all of the runs there was some slow oscillation. This 
observation will be commented upon in discussing the results. 

PROCEDURE 

The experimental run conditions were established, generally, by de- 
ciding on a void fraction and then varying the gas and liquid rates to 
give that void fraction over a broad range of total flow rates. In 
actual practice the data points were set by watching the indication on 
the momentum measurement and then adjusting the gas and liquid flows so 
that the void measurement remained constant. From this procedure, a con 
stant void curve was generated. This led eventually to a family of con- 
stant void curves covering a broad flow rate range. 

All of the data were recorded on strip charts. Average values for 
each measurement were estimated and these values were utilized in the 
computerized calculation of the data. . 

CALCUL AT I ONAL PROCEDURE 
One-Dimensional Model (Measured a) 

The momentum' of the two phase, two component flow discharging from 
the test section was computed using a one-dimensional, two-phase equa- 
tion. (ref. 1 or 2). The equation for the discharge force (rate change 
of momentum) is 

2 2 
- _ W w ' W N 

F CAL 8 c P w A( 1 - a) + g c P N Aa ■ ^ 

where a is an experimentally determined value in this case. The force 
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produced by the nitrogen contribution in equation (2) was negligible 
compared to the water contribution. Consequently, 

W 2 

— w 


CAL g p A(1 - a) 


( 3 ) 


c w 


or in terms of the liquid Reynolds number 


■F,, 


2 2 2 
RPy, 
w L 


CAL 16g p A(1 - a) 
°c w 


(A) 


The results from this calculation were compared with the experimentally 
measured impulse obtained from the momentum cage-load cell arrangement 
at the discharge end of the test section. 

Slip Ratio 

The ratio of the average vapor to liquid velocity in a two phase 
system is commonly called the slip ratio, k. This ratio is related to 
the quality and void fraction by the expression 

Pt 


, L x 1 - a 

k = — 

p 1 - x a 
v 


(5) 


where the quality x is defined as 


W- 


x = 


w + w 

w v 


( 6 ) 


Homogeneous Model (No Slip) 

As previously discussed, the simplest of the flow models is the 
homogeneous model in which both the liquid and vapor move at the same 
velocity (no slip). In this case, the homogeneous void fraction is 
given as 


HOMO , p 

1 

x p L 


( 7 ) 
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The calculated force for the homogeneous model is found by substituting . 
in equation (2) the values of a computed from equation (7). 

Anticipated Error in Force Estimates 
From equation (2) it is apparent that the force results will be very 
sensitive to the void fraction measurements. If the difference between 
the actual void fraction and the measured void fraction is 

a “ a MEAS = Aa 


then the percentage error in the value of F„ AT /F..„ . c would be 

UAL rlLAb 

Act 

percent error = — — x 100 (8) 

From the data shown in figure 3, Aa is estimated to be 0.05. Conse- 
quently, for values of a ^ 0.25 the computed force will have an uncer- 
tainty of at least ±6 percent. For a ^ 0.5, the uncertainty would in- 
crease to at least 10 percent. This uncertainty will be represented by 
vertical lines in the results, to be presented (figs. 7 and 8). 


DISCUSSION OF RESULTS 
a = 0 

For the all liquid case (a -'0), a comparison of the calculated to 
the measured force appears in figure 4 while the ratio of the calculated 
to measured force as a function of liquid velocity (Reynolds number) is 
shown in figure 5. As seen in figures 4 and 5, the calculated measured 
values are within 6 percent of each other. Some of this 6 percent devia- 
tion is experimental error and some could be attributed to changes in the 
velocity profiles over the range of operating conditions. The capacity 
of the load cell, 22.4 newtons (5 lb^) , limits the maximum flow and thus 
the upper value. of the Reynolds number of 100 000. The lowest Reynolds 
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number (40 000) is determined by the lower limit to which the pump can 
be throttled. 

Slip Ratio 

The slip ratios are shown in figure 6. Data were taken for nominal 
void fractions of 0.25 and 0.5. The computed slip results are listed in 
table I. Some of the scatter in the slip ratio results from the meas- 
ured variations in a. The void fraction varies from 0.45 to 0.55. In 
addition, it was estimated that the data have an uncertainty of about, 

±20 percent. 

As seen in figure 6, the slip ratio is nominally one when the void 
fraction is 0.25, but the slip ratio deviates significantly from one at 
the higher void fraction where a -0.5. 

Comparison of Variable Slip (Measured a) and Homogeneous Models 

For a :l 0.25, the ratio of the calculated to the measured force for 
both the homogeneous and variable slip (measured a) models are in very 
close agreement with the calibration line, as seen in figure 7. Thus, 
both of these one-dimensional two phase models appear to be applicable 
over the entire range of Reynolds numbers tested. This is not surpris- 
ing since the slip ratios for both models turned out to be near unity. 
When a 2 . 0,. 5, the slip ratios of the two models differed and they pre- 
dicted different force values. 

The Reynolds number in the abscissa of figure 7 is an all-liquid 
number. The data shown in figure 7 could just as well have been plotted 
against liquid velocity. 

For a ^.0.5, significant deviations between the one-dimensional 
variable slip, no slip (homogeneous), and the measured impulse force are 
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seen to occur in figure 8. In particular, a steep drop in the one- 
dimensional variable slip representation is seen in figure 8 near a 
Reynolds number of 0.4 X 10"\ For these low Reynolds numbers the variable 
slip calculation is in error. For larger Reynolds numbers, the variable 
slip representation seems to be approaching the calibration line, which -• 
indicates that this model is accurate at high liquid velocities (Reynolds 
numbers). This observation is in agreement with reference 5. The drop 
off in the variable slip model at high Reynolds numbers as suggested in 
reference 3 was not seen. However, such a trend might occur at higher 
void fraction. 

The homogeneous model fortuitously falls along the unity line in 
figure 8; which is in close agreement with the results of reference 3. 

The homogeneous model prediction is approximately 6 percent high, which 
for practical estimation purposes is quite satisfactory. It should be 
pointed out, however, that the predicted homogeneous values represent 
variable a values. The experimental data and the variable slip model 
are based on a constant a, which was measured by gamma attenuation. 

The a which one could calculate from the homogeneous expression would 
not. agree generally with the experimentally-determined value. 

One-Dimensional Variable Slip (Measured a) Correlation 

The drop off in the deviation between the computed variable slip 
force and the measured force at low Reynolds numbers. (~40 000) will now 
be discussed. 'Andeen and Griffith (ref. 3) present three possible 
reasons which might account for the above deviations at low Reynolds 
numbers: (1) - two-dimensional velocity profiles, (2) unsteady fluctua- 

tions Such as slugging, and (3) turbulent fluctuations. A fourth possi- 
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ble explanation (4) could be a transition in flow patterns (bubbly, 
annular, etc.)* In an actual situation, a combination of the above 
mechanism could account for the deviation between prediction and experi- 
ment. Each of these reasons will be discussed in detail because any or 
all of the above reasons can cause the real force to be greater than the 
one-dimensional variable slip value as calculated in equation (2). 

Spatial Velocity Distribution 

Two dimensional velocity profiles will lead to higher estimates of - 
the impulse force than would the simple one-dimensional models using an 
average channel velocity. 

For the case where a = 0 (all water) , if a power law velocity pro- 
file assumed for a single phase system is 

V “ V MAX (y/n)1/n (9) 

then- the ratio of the one-dimensional to the two-dimensional force is 





For the Reynolds number range considered herein, a value of n = 7 
is appropriate (ref. 7). In this case the above ratio is about 0.98. 

The two-dimensional impulse force is two percent greater than that cal- 
culated for the simple one-dimensional case. Physically, since the im- 
pulse force is proportional to the velocity squared, the fluid particles 
with velocities ■ greater than the average will produce a larger contribu- 
tion to the total force than those below the average. 

Could two-dimensional effects produce the high drop off shown in 
figure 8 for a = 0.5? For this to occur, the value of n in the two- 
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phase fully dispersed flow would of necessity be much smaller than that\ 
normally observed in turbulent flow of a single phase fluid. This is not 
expected. However, for other two-phase flow regimes, such as annular, 
significant changes in n are possible. 

Time Dependent Variations 

As mentioned earlier, unsteady fluctuations such as slugging or 
flow regime transitions will decrease the ratio of the one-dimensional 
calculated force to the measured force recorded by the load cell. The 
recorded output of the load cell showed the measured force to be peri- 
odic; consequently, time-averaged variations in the velocity must be 
considered. Figure 9 displays the amplitude of the load cell oscillation, 
F amp’ aS 3 ^ unct ^ on th e Reynolds number. As seen in figure 9, the 
force measurement exhibits a significant oscillatory behavior and the 
amplitude of the oscillations increase with increasing void fraction and 

<7 

decreasing Reynolds number. This is a reflection of a changing flow 
pattern. At the low void fraction a = 0.25, the flow configuration is 
generally one of bubbles in a continuous liquid stream - bubbly flow. ..As 
the void fraction increases, the bubbles coalesce into large gas volumes 
and the flow pattern might be better characterized as alternate slugs of 
gas and liquid-slug flow. With this observation, it is interesting to 
examine analytically what a large scale oscillation will do to the re- 
sults. 

The velocity oscillations which occur near the exit of the tube pro- 
duce the periodic output of the load cell. Unfortunately, the actual 
magnitude of the velocity oscillations are unknown. Conceivably, the 
amplitude of the velocity oscillation could be quite large in comparison 
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to the mean velocity v. 

Since the force exerted by the gas is negligible, we will consider 
only the liquid phase with cross sectional area A in the following 
calculations. Furthermore, we assume 

v = v(l + a cos tot) (11) 

where "a" times v represents the amplitude of the oscillation. The 
flow area will vary with the velocity. A detailed look at the mo- 

mentum and continuity equation will be required to determine the rela- 
tionship between v and A^. For the purpose of simplicity, however, 
we will assume A^ to be a constant. 

To preclude reverse flow, the amplitude factor a is assumed to 
have a limiting value of 1. During the course of the experiment no 
visual indications of flow reversal were seen. We will now investigate 
what effect the amplitude "a". will have on the calculated average force. 

The average force for one cycle (period = 2ir/m) can be calculated 
as follows. 



( 12 ) 


Substituting equation (11) into equation (12) and integrating gives 


t) < 13 > 

As seen in the above equation, the average force will be independent of 
frequency; however, the response of the measuring system will, of course, 
depend on. the frequency. Noting that 


_ PAjV 
F = L 


Sc 


1 + 


A^ = A(1 - a) 


(14) 
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and taking the ratio of equation (2) to equation ( 13 ) 


F 

CAL 1 



( 15 ) 


As previously mentioned, for the limiting case of a equal to 1, the 
lower limit to the ratio given by equation (15) would be 
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As seen. in figure 8, there is an indication that the ratio of the 
calculated to impulse force approaches the 2/3 value. Some detailed 
measurements of the velocity leaving the tube would be required to deter- 
mine if the 2/3 value suggested by equation (16) represents a lower limit 
to the data or if this is just a fortuitous occurrence. This calculation 
leaves us with the impression that flow oscillations could drastically 
affect the momentum results. Unfortunately ; we cannot draw any strong 
conclusions regarding this possibility as it applies to. the data reported 
herein. 


A word of caution is necessary at this point. Figure 10 displays a 
ratio, of the amplitude of the load cell oscillation, to the mean force 
amplitude, F, as a function of Reynolds number. As seen in figure 10, 
the ratio increases significantly at the low Reynolds numbers where the 
drop off occurred in figure 8. From the results in figure 10, we might 


expect that the ratio of F /F . c for . a = 0.25 in figure 7 should 

L-ALi MiiAo 


display a drop off similar to that shown in figure 8 for a = 0.5. This 


is contrary to the conclusion drawn from figure 9. We suspect, however. 


that because of damping and frequency response of the momentum cage, the 
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magnitude of oscillation for the a. = 0.5 data should be higher. Thus, 
it may be fortuitous that both sets of data fell on the single curve dis- 
played in figure 10. In our opinion, it is safer to draw conclusion from 
figure 9 rather than figure 10. 

i > 

Turbulent Fluctuations 

Turbulent fluctuations are somewhat similar to the time dependent 
fluctuations just discussed, if for example we assume 

v = v + v' (17) 

where v' is a fluctuation random velocity component. In general, v T 
could be represented by an infinite Fourier series; however, as rough 
approximation for turbulent effects, we might assume 

v = v(l + v' cos wt) (18) 

MAX 

where a crude representation of largest intensities seen in the 

channel. At any position in the channel, Vj]^ will generally be less 
than 0.15. Substituting this value into equation (15) indicates that the 
turbulent fluctuations may account for less than 1 percent of the drop 
observed in the calculated to measured force ratio. 

Flow Patterns 

The experimental flow regime in the test was nominally in the bubbly 
regime as indicated by the position of the experimental data on the Baker 
plot shown in figure 10. However, since in reality, the separation of 
bubbly from slug flow actually occurs over a broad range, a flow transi- 
tion might be beginning at the lower Reynolds numbers. 

CONCLUSIONS 

For the range of Reynolds numbers and flow regime investigated, the 
following conclusion can be made. 
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1. The one-dimensional variable slip (measured a) and the no-slip 
homogeneous predictions of the two phase momentum are valid for void 
fractions equal to or less than 25 percent. 

2. For void fractions of about 50 percent, the variable-slip model 
is valid for large Reynolds numbers. However, for Reynolds numbers below 
50 000 a significant deviation between experiment and prediction is seen. 
The homogeneous no-slip model seems to predict the impulse force for this 
higher void fraction over the entire Reynolds number range. This obser- 
vation, which is in agreement with the experimental results in refer- 
ence 3, is due to the over prediction of the one-dimensional momentum by 
the homogeneous model in which the slip is assumed to be unity. As 
illustrated by the data, the homogeneous model is a versatile and useful 
technique for one-dimensional design calculations in the range of 
Reynolds numbers investigated. 

LIST OF SYMBOLS 
cross sectional area of duct 
cross sectional flow area of liquid 
amplitude of velocity fluctuation 
average calculated force, equation (2) 
average force for one-dimensional flow (slug) 
average force for two-dimensional flow 
average measured force 
time averaged force 

2 

nitrogen flow rate, lb m/hr ft for Baker plot 

2 

water flow rate, lb m/hr ft for Baker plot 


A 

*L 

a 

F 


CAL 


ID 


2D 


MEAS 




w 


gravitational constant 
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h l/2-height of channel 

P perimeter 

R liquid Reynolds number 4W /Pu 

w w w 

n turbulent power law coefficient 

v velocity 

v average velocity 

v,.„ maximum velocity 

MAX 

v' fluctuating turbulent component of v 

v' maximum value of v' 

MAX 

nitrogen flow rate 

W v vapor flow rate 

W water flow rate 

w 

x quality, equation (6) 

y height coordinate 

y y/ h 


a 

Aot 

a HOMO 

a MEAS 

X 



void fraction 

01 " “meas 

homogeneous void fraction 
measured void fraction 

i/p 

equals [ (p^/0 . 075) /p^/62. 3 ] (for Baker plot) 

viscosity of water, centipoise for Baker plot 

density 

liquid density 

3 

density. of nitrogen, lb m/ft for Baker plot 
vapor density 

3 

density of water, lb m/ft for Baker plot 
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* 

+f 

*o 

rp 


1 . 


2 , 


3. 


4. 


5. 


6 . 


7. 


gamma signal 

gamma signal for all water 

empty gamma signal for all nitrogen 


equals (73/a) 
angular velocity 


y L (62.3/p L )‘ 


-1/3 


(For Baker plot) 
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Figure 1. -Momentum pressure drop apparatus. 


Figure 2. -Momentum cage and load cell. 
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Figure 3. - Comparison of theoretical and calibrated void 
as determined from gamma-ray attenuation experiment. 
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Figure 4. - Comparison of measured to calculated force 
for zero void fraction - all liquid. 
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Figure 5. - Ratio of force calculated from one di- 
mensional assumption to measured force for 
a = 0. 
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Figure 6. -Measured slip ratios. 
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gure 7. - Ratio of force calculated from one dimensional assumption 
to measure force for a = 0. 25. 
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Figure 8. - Ratio of force calculated from one dimensional assumptions to 
measured force for a - 0.5. 



Figure 9. - Amplitude of measured force for 
a = 0.25 and 0.5. 



Figure 10. - Ratio of amplitude to measured 
force. 



Figure 11. - Flow regimes of data as designed by baker 
plot. 
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